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ABSTRACT 


Improvised  Explosive  Device  (IED)  is  used  as  a  strategic  weapon  of  choice  and 
continues  to  be  a  threat  both  globally  and  domestically.  One  of  the  deadly  devices  in  this 
arsenal  is  the  Explosively  Formed  Projectile  (EPF).  This  study  develops  methodology  for 
modeling  and  simulation  of  armor  plates  to  survive  EFP  threats.  The  EFP  effects  are 
modeled  as  a  pressure  or  blast  wave  using  compressible  isentropic  conservation  equations 
to  get  pressure  loadings.  The  thermal  effects  are  modeled  as  temperature  intensities  and 
resulting  transient  heat  transfer  analysis  is  conducted  to  obtain  temperature  distribution. 
The  kinetic  loads  are  modeled  as  high  initial  velocities  applied  to  the  plate.  The  combined 
mechanical  and  thennal  loading  is  analyzed.  The  design  space  is  generated  for  varying 
materials  properties  and  thicknesses  as  parameters.  Laminated  composite  and  ortho  tropic 
composites  are  also  used  in  addition  to  special  high  strength  and  high  stiffness  generic 
alloys.  The  analysis  is  done  using  both  two-dimensional  plate  theories  as  well  as  three- 
dimensional  transient  dynamic  analysis.  The  results  are  presented  showing  maximum 
stresses  and  deformations  for  different  combinations  of  materials  and  thicknesses.  The 
results  also  indicate  the  need  to  use  three-dimensional  analysis  for  designing  survivable 
annor.  Some  recommendations  are  made  for  further  studies. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

A.  IMPROVISED  EXPLOSIVE  DEVICE  BACKGROUND . 1 

1.  Lessons  from  Iraq . 1 

B.  SURVIVABILITY . 2 

1.  Explosively  Formed  Projectile . 3 

II.  THEORY . 5 

A.  THEORETICAL  MODEL  OF  IED  EXPLOSION  AND  IMPACT . 5 

1.  Pressure  Wave . 5 

2.  Explosively  Formed  Projectile . 6 

B.  MATERIAL  ENERGY  ABSORPTION . 6 

1.  Deformation . 6 

C.  TRANSIENT  VIBRATION . 9 

D.  THERMAL  ANALYSIS . 10 

1.  Heat  Transfer . 10 

2.  Thermal  Stresses . 11 

III.  FINITE  ELEMENT  FORMULATION . 13 

IV.  PROCEDURE . 19 

V.  RESULTS . 21 

A.  KINETIC  LOADING  RESULTS . 2 1 

1.  Kinetic  Loading  Results  of  Various  Materials  at  10cm 

Thickness  in  2D . 21 

2.  Kinetic  Loading  Results  of  Steel  at  Varying  Thickness  in  2D . 29 

B.  COMPOSITES  WITH  INTER-LAMINAR  SHEAR . 36 

C.  THERMAL  LOADING  RESULTS . 37 

D.  3-D  THERMAL  LOADING  RESULTS . 41 

E.  3-D  COMBINED  LOADING  OF  ALUMINUM . 45 

VI.  CONCLUSION . 51 

APPENDIX . 53 

A.  ARMOR  PLATE  CONSTRUCTION . 53 

1.  Our  First  Step  is  to  Create  a  New  Data  Base . 53 

2.  Next  Set  the  Analysis  Preference . 53 

3.  The  Geometry  of  the  Beam  will  Beam  Determined  Next . 53 

4.  Specify  the  Boundary  Conditions  Next . 54 

5.  The  Finite  Element  Mesh  is  Specified  Next . 55 

B.  ARMOR  PLATE  LOAD  CASE  CONSTRUCTION . 56 

1.  The  Load  Case  is  Created  Next . 56 

2.  The  Loads  are  Specified  Next . 56 

C.  CREATING  ARMOR  MATERIAL . 58 

1.  The  Materials  are  Specified  Next . 58 

vii 


2.  The  Properties  are  Assigned  Next . 60 

D.  KINEMATIC  LOAD  ANALYSIS . 61 

1.  Analysis  is  Performed  Next . 61 

2.  Accesses  the  Results . 63 

E.  THERMAL  LOAD  ANALYSIS . 63 

F.  THERMAL  STRESS  ANALYSIS . 63 

1.  Accesses  Thermal  the  Results . 63 

2.  Create  a  Spatial  FEM  Based  on  the  Temperature  Profile . 64 

3.  Change  the  Analysis  Type . 64 

4.  Define  a  Temperature  Load . 64 

LIST  OF  REFERENCES . 67 

INITIAL  DISTRIBUTION  LIST . 69 


viii 


LIST  OF  FIGURES 


Figure  1 .  HUMV  destroyed  by  IED . 1 

Figure  2.  Diagram  of  simple  IED  attack . 2 

Figure  3.  Diagram  of  EFP  detonation  and  projectile  shaping . 3 

Figure  4.  Displacement  under  uniaxial  stress . 7 

Figure  5.  Work  done  by  uniaxial  stress  [2] . 8 

Figure  6.  Deformation  due  to  pure  shear  [2] . 9 

Figure  7.  Stress  distribution  across  an  aluminum  plate . 22 

Figure  8.  Stress  distribution  across  a  steel  plate . 23 

Figure  9.  Stress  distribution  across  a  high-modulus  graphite  epoxy  plate, 

homogeneous . 23 

Figure  10.  Stress  distribution  across  a  high-strength  graphite  epoxy  plate, 

homogeneous . 24 

Figure  11.  Stress  distribution  across  a  multi-material  plate,  homogeneous . 24 

Figure  12.  Plot  of  Von  Mises  Stresses  of  tested  materials . 25 

Figure  13.  Deflections  across  an  aluminum  plate . 26 

Figure  14.  Deflections  across  a  steel  plate . 26 

Figure  15.  Deflections  across  a  high-modulus  graphite  epoxy  plate,  homogeneous . 27 

Figure  16.  Deflections  across  a  high-strength  graphite  epoxy  plate,  homogeneous . 27 

Figure  17.  Deflections  across  a  multi-material  plate . 28 

Figure  18.  Plot  of  deflections  of  tested  materials . 29 

Figure  19.  Stress  distribution  across  a  4  cm  thick  steel  plate . 30 

Figure  20.  Stress  distribution  across  a  6  cm  thick  steel  plate . 30 

Figure  2 1 .  Stress  distribution  across  an  8  cm  thick  steel  plate . 31 

Figure  22.  Stress  distribution  across  a  10  cm  thick  steel  plate . 3 1 

Figure  23.  Stress  distribution  across  a  12  cm  thick  steel  plate . 32 

Figure  24.  Displacements  across  a  4  cm  thick  steel  plate . 32 

Figure  25.  Displacements  across  a  6  cm  thick  steel  plate . 33 

Figure  26.  Displacements  across  an  8  cm  thick  steel  plate . 33 

Figure  27.  Displacements  across  a  10  cm  thick  steel  plate . 34 

Figure  28.  Displacements  across  a  12  cm  thick  steel  plate . 34 

Figure  29.  Chart  of  steel’s  performance  at  different  thicknesses . 35 

Figure  30.  Tsai-Wu  criteria  stresses  in  high-modulus  graphite . 36 

Figure  3 1 .  Tsai-Wu  criteria  strains  in  high-modulus  graphite . 37 

Figure  32.  Temperature  distribution  on  an  aluminum  plate . 38 

Figure  33.  Temperature  distribution  on  a  steel  plate . 38 

Figure  34.  Von  Mises  stress  distribution  caused  by  initial  temperature  on  aluminum . 39 

Figure  35.  Von  Mises  stress  distribution  caused  by  initial  temperature  on  steel . 39 

Figure  36.  Displacements  from  initial  temperature  on  aluminum . 40 

Figure  37.  Displacements  from  initial  temperature  on  steel . 40 

Figure  38.  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  0.0005 . 41 

Figure  39.  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  0.507 . 42 

Figure  40.  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  1 .003 . 42 


IX 


Figure  4 1 .  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  1 .499 . 43 

Figure  42.  Heat  flux  acting  on  a  3d  aluminum  plate  at  t  =  0.0005 . 43 

Figure  43.  Heat  flux  acting  on  a  3d  aluminum  plate  at  t  =  1 .499 . 44 

Figure  44.  Temperature  gradient  across  a  3d  aluminum  plate  at  t  =  0.0005 . 44 

Figure  45.  Temperature  gradient  across  a  3d  aluminum  plate  at  t  =  1.499 . 45 

Figure  46.  Temperature  gradient  across  a  3d  aluminum  plate  at  t  =  0.5  sec . 46 

Figure  47.  Maximum  von  Mises  stress  on  a  3d  aluminum  plate  at  t  =  0.5  sec . 46 

Figure  48.  Maximum  displacements  on  a  3d  aluminum  plate  at  t  =  0.5  sec . 47 

Figure  49.  Maximum  von  Mises  stress  on  a  3d  steel  plate  at  t  =  0.5  sec . 48 

Figure  50.  Maximum  displacements  on  a  3d  steel  plate  at  t  =  0.5  sec . 49 


x 


LIST  OF  TABLES 


Table  1 .  Summary  of  key  data  for  steel  of  various  thickness  x  3  m  x  2  m . 51 

Table  2.  Weight  of  3  m  x  2  m  x  0. 1  m  armor  plate  made  of  listed  materials . 52 

Table  3.  Properties  of  various  tested  materials . 59 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


ACKNOWLEDGMENTS 


I  would  like  to  thank  Professor  Ramesh  Kolar  for  all  the  help  he  gave  me.  His 
extensive  knowledge  made  it  possible  to  quickly  troubleshoot  any  model  problems, 
making  it  possible  to  complete  the  present  work  on  time. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xiv 


I.  INTRODUCTION 


A.  IMPROVISED  EXPLOSIVE  DEVICE  BACKGROUND 

1.  Lessons  from  Iraq 

Improvised  Explosive  Devices  (IEDs)  are  the  single  greatest  cause  of  causalities 
and  destruction  of  equipment  in  Iraq,  Afghanistan,  and  elsewhere  as  shown  in  Figure  1. 
IEDs  can  be  detonated  in  a  multitude  of  ways  from  command  detonation  with  a  cell 
phone  to  a  command  detonation  from  a  sensor  to  a  simple  trip  wire.  The  resulting 
explosion  produces  a  pressure  wave  that  causes  sever  structural  damage  to  a  vehicle  and 
can  in  some  cases  lift  the  vehicle  off  the  ground  and  flip  it  over.  More  advanced  IEDs 
also  use  an  explosively  formed  projectile  (EFP)  to  penetrate  vehicle  armor.  The 
explosively  formed  projectile  impacts  the  target  with  a  molten  mass  of  metal  at  velocities 
of  1 .5  to  2  km  per  a  second. 


Figure  1.  HUMV  destroyed  by  IED 
1 


One  solution  to  defend  against  IEDs  and  EFPs  is  military  vehicles  must  be 
designed  to  mitigate  the  effects  of  EFPs  to  improve  survivability  of  the  vehicle  by  the  use 
of  advanced  forms  of  annor.  The  old  fashioned  method  of  just  adding  inches  of  steel  are 
not  practical  and  would  result  in  viability  and  maintenance  issues  with  the  vehicles. 
Additionally  such  a  heavy  vehicle  would  have  a  difficult  time  patrolling  an  urban  area 
with  its  narrow  streets  and  restrictive  weight  limits  on  the  roads  and  bridges. 

B.  SURVIVABILITY 

Survivability  is  simply  the  ability  to  of  a  military  vehicle  to  go  out  on  a  mission 
and  come  back  to  do  it  again.  Survivability  is  composed  primarily  of  two  elements, 
susceptibility  and  vulnerability.  Susceptibility  is  the  inability  of  a  vehicle  to  avoid  being 
damaged  by  the  hostile  environment.  Vulnerability  is  the  inability  of  a  vehicle  to 
withstand  damage  caused  by  the  hostile  environment.  This  paper  will  focus  on  mitigating 
the  effects  of  an  EFP  by  reducing  the  vulnerability  of  military  vehicles  by  modeling  a 
typical  IED  explosion  and  developing  methodology  for  modeling  and  simulation  of 
armor  that  will  withstand  that  attack  allowing  the  crew  to  evacuate  and  fight  on.  A 
simple  IED  attack  looks  like  Figure  2. 


The  Basie  “No  Frilfs"  IED  Attack 


Mil/Civ  Convoy 


Direction 
of  Traffic 


IED  Placed  on  Shoulder 

\\ 

3 

Median  Strip 

Median  Strip 

Median  Strip 

l=> 


Typical  Iraoi  4  lane  highway 


Varia  tion:  IED  Placed  in  the  Median  Strip 
Figure  2.  Diagram  of  simple  IED  attack. 


2 


1.  Explosively  Formed  Projectile 


An  Explosively  Formed  Penetrator  (EFP)  is  a  class  of  IEDs  with  a  shaped 
warhead  designed  to  penetrate  armor  effectively  at  stand-off  distances.  An  EFP  system 
comprises  three  major  components:  1)  an  arming  device,  2)  a  trigger  device  and  3)  the 
explosive  device.  Figure  3  shows  a  diagram  for  an  EFP.  The  explosive  device 
component  consists  of  four  elements:  a  concave  metal  liner  (warhead),  case,  base  plate, 
and  explosive  charge  (propellant).  The  case  holding  the  explosive  is  generally 
cylindrical,  fabricated  from  commonly  available  materials  (e.g.  PVC  pipe,  steel  pipe), 
with  the  forward  end  closed  by  a  concave  copper  or  steel  disk-shaped  liner  to  create  a 
penetrator.  Generally  military  grade  explosive  is  loaded  behind  the  metal  liner  to  fill  the 
casing.  A  blasting  cap  initiator  is  placed  through  a  hole  centered  in  the  base  plate.  Upon 
detonation,  the  explosive  projects  the  liner  to  form  a  projectile  in  a  direction  in  line  with 
the  EFP  at  a  speed  well  over  one  kilometer  per  a  second,  depending  on  the  design  and  the 
type  of  explosive  used.  Detonation  is  controlled  by  command  wire,  radio  control  (RC), 
or  remote  arming  with  trigger.  EFPs  can  be  deployed  singly,  in  pairs,  or  in  arrays,  and 
camouflaged  in  painted  foam  to  look  similar  to  natural  rocks  or  with  man  made  artifacts 
(e.g.  trash)  depending  on  the  tactical  situation. 


Figure  3.  Diagram  of  EFP  detonation  and  projectile  shaping. 
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In  order  for  the  crew  of  a  vehicle  to  survive  and  evacuate,  the  vehicle’s  armor 
would  have  to  with  stand  the  resulting  pressure  wave  and  prevent  the  EFP  from  entering 
the  crew  cabin. 
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II.  THEORY 


A. 


THEORETICAL  MODEL  OF  IED  EXPLOSION  AND  IMPACT 


1. 


Pressure  Wave 


In  this  section  a  simple  IED  explosion  will  be  examined  from  the  point  of 
explosion  to  the  point  of  impact,  on  a  rectangular  plate  of  armor.  When  an  IED  is 
detonated;  the  resulting  explosion  produces  a  pressure  wave.  Typically  the  explosion  is 
modeled  [1]  as  a  compressible,  time-dependent,  conservation  equations  for  the  mass 
density,  momentum  density  and  energy  density.  These  equations,  written  in  their 
Eulerian  form  and  in  one-dimensional  spherical  coordinates  are  [1]: 

1  dr2  pv 


dp 

dt 


dpv 

dt 


r  dr 

2 


de 

dt 


1  dr  pw  dp 
r2  dr  dr 
1  dr2ev  1  dr2vp 


(1) 

(2) 

(3) 


r  dr  r  dr 

where  p  is  the  density,  v  the  fluid  velocity,  p  the  pressure  and  e  the  volumetric  density 
of  the  total  energy,  v  and  t  are  the  radial  coordinate  and  time  respectively.  The  total 
energy  density  e  is  defined  as  the  sum  of  the  internal  energy  and  kinetic  energy.  Using 
equation  4  the  work  done  by  the  gas  during  an  isentropic  expansion  can  be  expressed  as 


wis=pyx 


1 

f  Wi-'fn  1 

1- 

Po_ 

/l-1 

[PiJ 

(4) 


where  P0  is  the  ambient  temperature  and  P}  is  the  expected  over  pressure. 

The  peak  pressure  ps  and  distance  R  are  non-dimens ionalized  using  Sachs’s 
scaling  relationship s.  [  1  ] 

Ps  ={Ps~  Po)!  Po  (5) 
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R  = 


R 


/  \  1/3 

'W„  ' 


VPoJ 

The  peak  over  pressure  can  then  be  calculated  using  equation  6. 

Ps  ~Po  =  P  Po 


(6) 


(7) 


This  resulting  overpressure  is  applied  as  a  pressure  force  on  all  subsequent 
simulations  of  the  IED  effects  on  the  annor  plate.. 


2.  Explosively  Formed  Projectile 

In  this  study,  the  EFP  effects  are  modeled  in  two  steps.  The  thennal  effects  are 
modeled  as  high  temperature  flux  on  the  plate.  The  kinetic  energy  is  modeled  as  an 
applied  initial  velocity  to  the  plate.  The  equations  of  motion  may  be  expresses  as  follows 

[m]{;}+[k]{*}={/>)  <8) 

where  [M]  is  the  mass  matrix,  [K]  the  stiffness  matrix,  {x}  the  displacement  vector,  and 
{P}  the  loading  vector.  Direct  numerical  integration  method  is  used  in  computing  the 
response  to  the  combined  loads.  [8] 

B.  MATERIAL  ENERGY  ABSORPTION 


1.  Deformation 

Materials  absorb  mechanical  energy  through  the  process  of  deformation.  At  first 
the  deformation  is  elastic  and  any  changes  in  shape  will  revert  to  the  original  shape  once 
the  load  is  removed.  After  elastic  deformation  comes  inelastic;  any  changes  in  shape  will 
be  pennanent.  The  amount  of  deformation,  both  elastic  and  inelastic,  that  a  material  can 
undergo  is  called  ductility.  The  more  ductile  the  material  the  more  it  can  deform  and  the 
more  energy  it  will  absorb. 
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The  work  done  by  external  forces  in  causing  deformation  is  stored  within  the 
body  in  the  form  of  strain  energy.  In  an  ideal  elastic  process,  no  dissipation  of  energy 
takes  place  and  all  the  stored  energy  is  recoverable  upon  unloading. 

The  temperature  distribution  is  used  to  conduct  a  transient  heat  transfer  analysis 
resulting  in  temperature  distributions  in  the  plate  as  a  function  of  time.  In  the  next  step, 
these  temperatures  are  applied  as  thermal  loads  together  with  blast  loads  and  initial 
velocity  to  perform  a  transient  dynamic  analysis  resulting  in  displacements,  strains,  and 
stresses  caused  in  the  armor  plate. 


Figure  4.  Displacement  under  uniaxial  stress  from  [2], 


The  net  work  done  on  an  element  can  be  calculated  by  [2]: 

f p,,  A 


dW  =  dU  =  [  <Jxd  —  dx  dydz  =  [  <Jxdsx  ( dxdydz )  = 

o  vSx  )  ‘  ‘ 


(9) 


o  V  /  o 

Note  that  dW  is  the  work  done  on  dx  dy  dz,  and  dU  is  the  corresponding  increase 
in  strain  energy.  Designating  the  strain  energy  per  a  unit  volume  as  Uo,  for  a  linearly 
elastic  material  in  one  dimension,  we  have 


U0  =  J  cjxdsx  =J  Esxdsx 


(10) 
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(11) 


After  integration,  Equation  10  yields: 


T  T  1  77  2  1  2  1  2 

Ur,  =  —Esx  =  -<jx£x  = - CT 

0  2  2  2E  * 


This  quantity  represents  the  shaded  area  in  Figure  5.  The  units  of  Uo,  strain 

a 

energy  density,  in  SI  are  joules  per  cubic  meter  (J/m  ),  or  pascals  (Pa). 


Figure  5.  Work  done  by  uniaxial  stress  from  [2], 

The  shaded  area  under  the  stress-strain  curve  represents  the  measure  of  a 
material’s  ability  to  absorb  energy  without  permanent  deformation,  and  is  defined  as  the 
modulus  of  resilience.  The  area  under  the  entire  stress-strain  curve  provides  a  measure  of 
a  material’s  ability  to  absorb  energy  up  to  the  point  of  fracture,  and  is  defined  as  the 
modulus  of  toughness. 

Shear  deformation  and  the  elastic  strain  energy  associated  with  it  is  now  analyzed 
by  considering  an  element  of  thickness  dz  subject  to  only  shearing  stress  ixy,  see  Figure  6. 
Note  that  the  shear  force  xxy  dxdz  causes  a  displacement  of  yxy  dy.  The  strain  energy  due 

to  shear  is  ( z v  dxdz  )  ( yxvdyj ,  where  the  factor  of  Vi  arises  because  the  stress  varies 

linearly  with  the  strain. 
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/  dx 

zr 


Figure  6.  Deformation  due  to  pure  shear  from  [2]. 

The  strain  energy  is. 

U0=-t  Y  =  —  r2  =—Gv2  (12) 

Because  the  work  done  by  xxy  accompanying  perpendicular  strains  yyz  and  yxz  is 
zero,  the  total  strain  energy  density  done  by  shear  alone  is  found  by  superposition  of  three 
tenns  identical  in  the  form  of  Equation  11. 

u0  +  W  +  Tx*rxz)  (13) 

C.  TRANSIENT  VIBRATION 

The  force  of  the  explosion  from  the  IED  will  impact  on  the  annor  plate  for  very 
short  finite  time  with  a  force  of  large  magnitude.  This  impact  of  force  application  is  best 
modeled  as  an  Impulse  Excitation  [3], 

F  =  jF(t)dt  (14) 

Because  F  dt  =  m  dv,  the  impulse  acting  on  the  armor  will  result  in  sudden  change 
in  the  armor’ s  velocity  equal  to  F  trn  without  appreciable  change  in  its  displacement. 
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Under  free  vibration,  the  undamped  single  degree  of  freedom  spring-mass  system  with 
initial  displacement  of  x  ( 0 )  and  initial  velocity  of  x  ( 0 )  given  by  the  equation 


x 


(<) 


(0)  . 

—  sm  cot 


co„ 


+  x(0)cos<y;,t 


(15) 


Hence,  the  response  of  a  spring  mass  system  initially  at  rest  and  excited  by  an  impulse  F 
is 


Jr  a. 

x  = - sin  cont  =  Fhit ) 

mcon 

where 


hit )  = - sin  cont 

mcon 


is  the  response  to  a  unit  impulse. 


(16) 


(17) 


D. 


THERMAL  ANALYSIS 


1. 


Heat  Transfer 


The  impact  of  the  copper  projectile  from  an  EFP  will  result  in  thermal  loads  to  the 
annor  plate.  For  the  purpose  of  this  research  the  thermal  effects  of  an  EFP  will  be 
modeled  as  an  initial  temperature  of  1400K  acting  at  the  same  point  as  the  initial  velocity. 
The  armor  design  requirements  need  to  dissipate  this  heat  energy  effectively  and 
efficiently.  We  assume  that  heat  will  transfer  throughout  the  armor  by  conduction  in 
accordance  with  Fourier’s  law  [4] 


q  =  -kV  T  =  -k 


dT  dT  ,  dT 

l - +7 - +  k - 

dx  dv  dz 


(18) 


dy  (7Z  J 

where  q  is  the  heat  flux,  k  is  the  thermal  conductivity,  (k  is  negative  to  indicate  the  flow 
of  heat  from  hot  to  cold)  T  is  the  temperature,  and  V  is  the  gradient  operator 
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2. 


Thermal  Stresses 


The  resulting  increase  in  temperature  throughout  the  annor  as  a  result  of  the 
temperature  flux  will  cause  the  armor  to  expand,  via  thermal  expansion 

st  =  aAT  (19) 

where  et  is  the  thennal  strain,  a  is  the  coefficient  of  linear  thennal  expansion,  and  AT 


is  the  change  in  temperature.  The  total  strains  on  the  armor  plate  in  two  dimensions 
would  then  become  [2] 


£x=-^(<rx-v<ry)  +  cAT 

£y  =  (o-y  ~  0(Jx )  +  aAT  (20) 


Notice  that  the  thermal  expansion  doesn’t  result  in  shear  distortion  for  isotropic  materials, 
where  as  in  a  composite  material  this  will  not  hold  true. 


Since  the  armor  is  constrained  on  all  its  edges  the  resulting  stress  in  terms  of 
strain  components  become  [2] 


a*  =  YIs(S-+USr) 

!J-  =  YZ^(£'+ve‘) 
t  =  Gy 

xy  /  xy 


EaAT 
1  -v 
EaAT 
1  -v 


(21) 


In  equation  (21),  the  first  term  corresponds  to  the  mechanical  loads  (blast,  initial 
velocity)  and  the  second  term  corresponds  to  the  thermal  effects. 
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III.  FINITE  ELEMENT  FORMULATION 


Finite  elements  represent  spatial  descretization  of  a  continuum.  As  such, 
however,  they  do  not  immediately  impose  nonlinearity.  When  nonlinearity  has  to  be 
taken  into  account  for  larger  displacements  and  /or  stress,  a  numerical  model  poses  new 
dimensions  to  the  descretization  in  addition  to  the  spatial  descretization.  That  is,  the 
descretization  is  applied  to  time,  load,  and  material  properties  by  using  piecewise  linear 
curves.  While  descretization  allows  approximate  solutions  by  numerical  methods,  it 
introduces  numerous  mathematical  singularities  which  may  complicate  computational 
processes.  Fortunately,  the  efficiency  of  modern  digital  computers  makes  it  feasible  to 
apply  complicated  computational  procedures  to  the  complex  systems  of  engineering 
problems. 


The  variational  principle  renders  the  system  governing  equilibrium  equations 
when  applied  to  a  functional  n  ,  representing  a  total  potential  of  a  continuum 

n  =  U-W  (22) 

where  U  is  the  strain  energy  of  the  system  and  W  is  the  potential  energy  of  the  external 
loads.  The  equilibrium  equations  can  be  obtained  by  invoking  the  principle  of  virtual 
work  or  the  Ritz  method 


8  FI  =  0  or 


(23) 


which  implies  that  the  total  potential  of  the  system  must  be  stationary  with  respect  to  the 
state  variables  (displacement)  for  equilibrium  to  be  ensured.  The  functional  IT  is  so 
called  because  it  involves  the  integral  of  implicit  functions  of  the  state  variables,  {u} . 


Considering  a  three-dimensional  continuum  for  a  nonlinear  problem,  the 
stationarity  condition  results  in 

J  <jj8  Eij  dV  =  | bf8  Si  dV  +  J  t.8 Si  dS  +  ^  pt8 uj  (24) 

V  VS  1 

Where  the  dots  and  8  denote  infinitesimal  increments  and  arbitrary  variations, 
respectively. 
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In  this  section  a  brief  overview  of  the  finite  element  approach  for  modeling 
transient  thermo-mechanical  response  is  presented  [8]. 

The  finite  element  method  can  be  characterized  by  the  following  features 
distinguished  from  the  conventional  Ritz  methods  or  the  matrix  method  for  frame 
structures: 

•  The  whole  region  of  the  system  is  divided  into  numerous  subdomains,  called 
finite  elements,  which  have  simple  geometrical  shapes. 

•  The  variational  process  is  limited  to  each  finite  element,  which  aggregates  into  a 
whole  region  when  assembled. 

•  The  admissible  displacement  field  within  each  element,  jwj,  can  be  expressed  in 

tenns  of  nodal  displacements  using  interpolation  functions  known  as  shape 
functions,  N, 

{«}  =  [#]{«}  (25) 

where  {w}  is  a  displacement  vector  consisting  of  all  nodal  points  of  the  element. 

The  strain  displacement  relations  for  the  element  can  then  be  established  in  terms  of 
nodal  displacements  using  the  shape  functions  in  Equation  26 

{^}  =  [£]{«}  (26) 

where 


and  the  element  matrix  [B]  consists  of  derivatives  of  the  shape  functions,  evaluated  at  the 
current  deformed  geometry.  MSC/NASTRAN  employs  an  approximated  Lagrangian 
approach  for  geometric  nonlinear  problems,  by  which  linear  strains  are  computed  in  the 
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updated  element  coordinate  system  in  order  to  eliminate  the  effects  of  the  rigid  body 
rotation  but  the  equilibrium  is  established  at  the  final  position  in  the  stationary  coordinate 
system. 

Equilibrium  equations  for  an  element  may  be  obtained  by  reducing  Equation  24 
after  the  substitution  of  Equations  25  and  26,  based  on  the  small  deformation  theory. 
Then  the  element  boundary  stresses  are  equivalent  to  the  nodal  forces  which  balance  the 
applied  external  loads 


(28) 

{F}'=J[S]»rfF 

(29) 

v 


and 

{pY  =\[Ni  {b)dV  +  \[N^  {t}dS  +  {p}  (30) 

V  s 

where  [/Vs  ]  is  an  appropriate  interpolation  function  for  the  traction  force.  Notice  that  the 
equilibrium  equation  for  an  incremental  load  may  be  expressed  as 

{>}  =  J[2?f  jaJjE  =  j>J  (31) 

where  jcrj  should  be  components  of  co-rotational  stress  which  is  independent  of  a  rigid 
body  rotation. 

The  element  stiffness  matrix  can  be  obtained  by  substituting  the  constitutive 
relations  into  Equation  32, 

{;}=[z>]{;}  P2) 

where 

[°}T  =  {^y<rzrxyTyzTa  )  (33) 
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and  [D]  is  a  tangent  material  matrix.  The  nodal  forces  of  an  element  can  then  be 
expressed  as 

|fJ  =pf  ja}r/E  =  [7f]ejwJ  (34) 

where  the  element  stiffness  is 

[K]  =  \[B^[D][B]dV  (35) 

V 

Notice  that  this  expression  represents  an  element  stiffness  due  to  material  stiffness 
without  geometric  nonlinear  effects. 

The  equilibrium  must  be  satisfied  in  the  whole  region  throughout  the  complete 
history  of  load  application.  Equilibrium  equations  for  the  global  discrete  system  are 
obtained  when  all  the  elements  are  assembled, 

(36) 

M  v  M 

where  Z  over  M  denotes  a  summation  over  all  elements.  For  the  incremental  process,  the 
equilibrium  equations  may  be  rewritten  as 

Zj[sfiCT-CT°}rfF= M  <37> 

M  v 

with 

M=ZM'-Zlb]VK  f3*) 

M  My 

where  jcr°j  represents  an  initial  stress  of  the  stress  state  at  the  preceding  load  step. 

Because  of  the  approximations  involved  in  the  interpolation  functions,  the  finite 
element  model  provides  an  approximate  solution  even  if  the  equilibrium  Equation  37  is 
satisfied  exactly.  Consequently,  the  differential  equations  of  equilibrium  are  not  satisfied 
exactly  even  for  linear  problems,  but  the  error  decreases  as  the  finite  element  mesh  is 
refined.  This  convergence  condition  is  required  and  ensured  by  element  formulations 
with  regard  to  the  element  convergence  criteria.  The  convergence,  however,  may  not  be 
monotonic  due  to  non-conforming  elements  or  reduced  integration.  By  virtue  of  the  non- 
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conforming  elements  and  the  reduced  integration,  the  finite  element  model  will  have 
added  flexibility  and  compensate  the  stiffening  effects  by  the  displacement  method.  With 
a  displacement  approach,  the  finite  element  model  is  generally  known  to  produce  a  stiffer 
structure  than  in  reality. 


17 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


18 


IV.  PROCEDURE 


In  the  present  study  the  engineering  analysis  tools  by  MSC  software  consisting  of 
PATRAN  and  NASTRAN  is  used  for  modeling  and  simulation.  MSC/PATRAN  is 
utilized  to  model  the  geometry  of  the  annor  plate,  the  materials,  thickness  properties, 
temperature  fluxes,  pressure  loads,  initial  conditions  and  boundary  conditions.  Here  in 
procedures  and  parameters  are  also  selected  in  this  module. 

The  input  deck  for  MSC/NASTRAN  is  generated  and  edited  if  necessary.  The 
MSC/NASTRAN  solver  is  used  to  obtain  the  solution  from  the  simulation.  The 
MSC/PATRAN  is  invoked  again  to  post-process  the  results  and  analyze  the  behavior  of 
the  annor  design. 

In  the  remains  of  this  chapter,  detailed  procedural  steps  to  model  and  simulate  the 
transient  thermal  and  transient  thermo-elastic  analysis  is  presented. 

The  armor  plate  is  modeled  as  a  3m  x  2m  plate  with  all  edges  clamped.  The 
temperatures,  initial  velocities,  pressure  loading,  material  properties,  thickness  properties, 
composite  laminate  properties  and  lay-up  are  used  as  design  parameters. 

The  armor  plate  was  modeled  as  follows. 

•  2-D  plate.  In  this  case  the  annor  plate  is  modeled  as  a  plate/shell  structure 
using  QUAD4  elements  of  the  analyzing  procedure. 

•  2-D  laminated  plate.  In  this  case  the  armor  plate  is  modeled  as  a 
laminated  composite  plate  with  each  laminate  modeled  as  an  orthotropic 
layer. 

•  3-D  plate.  In  this  case  the  armor  plate  is  modeled  as  a  three  dimensional 
structure  using  higher  order  tetrahedral  elements  (tetlO)  of  the  analysis 
procedure. 
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It  may  be  noted  that  for  laminated  composite  plate  analysis  only  blast  loads  and 
kinetic  impact  loads  were  considered  without  the  thennal  loads  due  to  software 
limitations.  Please  see  the  Appendix  for  detailed  procedures. 
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V.  RESULTS 


In  this  chapter,  the  results  for  thermal  analysis,  thermo-elastic  analysis  under 
thennal,  blast  and  kinetic  loads  are  presented.  It  is  noted  that  transient  heat  conduction 
analysis  is  perfonned  to  obtain  temperature  as  a  function  of  time.  These  temperatures  are 
then  applied  to  the  finite  elements  as  thennal  loads.  The  mechanical  loads  are 
superimposed  with  thennal  loads.  The  kinetic  loads  are  incorporated  as  prescribed  initial 
velocities.  A  transient  analysis  using  direct  solution  is  used  to  obtain  the  resulting 
dynamic  response  of  the  armor  plate.  Initial  linear  free  vibration  analysis  is  done  to 
select  appropriate  time  step  for  direct  solution. 

A.  KINETIC  LOADING  RESULTS 

1.  Kinetic  Loading  Results  of  Various  Materials  at  10cm  Thickness  in 
2D 

The  kinetic  loads  applied  to  the  annor  plate  include  the  pressure  wave  evenly 
distributed  across  the  top  surface  of  the  plate,  and  the  initial  velocity  applied  to  the  center 
of  the  plate.  The  results  are  presented  in  the  following  figures,  where  the  the  stresses  are 
in  Pascals  (Pa)  and  deflections  are  in  meters  (m). 
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Figure  7.  Stress  distribution  across  an  aluminum  plate. 
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Figure  8.  Stress  distribution  across  a  steel  plate. 
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Figure  9.  Stress  distribution  across  a  high-modulus  graphite  epoxy  plate,  homogeneous. 
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Figure  10.  Stress  distribution  across  a  high-strength  graphite  epoxy  plate,  homogeneous. 
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Figure  1 1 .  Stress  distribution  across  a  multi-material  plate,  homogeneous. 
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Max  Von  Mises  Stress 


Figure  12.  Plot  of  Von  Mises  Stresses  of  tested  materials. 

The  results  of  various  simulations  are  summarized  in  Figure  12.  The  x-axis  is 
used  to  identify  different  materials  and  y-axis  indicates  Von-Mises  stresses.  It  can  be 
seen  that  the  ASTM-A242  steel  had  the  highest  stress  while  the  graphite  composites  had 
the  lowest  with  the  high-strength  graphite  being  one  order  of  magnitude  lower  than  the 
rest. 
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Figure  13.  Deflections  across  an  aluminum  plate. 
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Figure  14.  Deflections  across  a  steel  plate. 
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Figure  15.  Deflections  across  a  high-modulus  graphite  epoxy  plate,  homogeneous. 
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Figure  16.  Deflections  across  a  high-strength  graphite  epoxy  plate,  homogeneous. 
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Figure  17.  Deflections  across  a  multi-material  plate. 
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Deflection 


Figure  18.  Plot  of  deflections  of  tested  materials. 

The  deflections(max)  for  avious  materials  is  summarized  in  Figure  18.  It  is  seen 
the  stresses  using  steel  has  the  greatest  deflections  while  the  graphite  composites  have  the 
least  and  the  high-strength  composite  is  an  order  of  magnitude  lower.  The  composites  in 
this  case  were  analyzed  as  orthotropic  material  in  NASTRAN. 

2.  Kinetic  Loading  Results  of  Steel  at  Varying  Thickness  in  2D 

The  results  for  the  armor  plate  under  kinetic  loading  are  presented  here.  The 
maximum  Von  Mises  stress  and  maximum  deflection  of  steel  as  varying  thickness  was 
found  to  help  determine  the  optimum  thickness  for  the  armor  plate.  The  kinetic  loads 
didn’t  change  only  the  thickness  was  changed  in  the  model. 
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Figure  19.  Stress  distribution  across  a  4  cm  thick  steel  plate. 
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Figure  20.  Stress  distribution  across  a  6  cm  thick  steel  plate. 
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Figure  2 1 .  Stress  distribution  across  an  8  cm  thick  steel  plate. 
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Figure  22.  Stress  distribution  across  a  10  cm  thick  steel  plate. 
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Figure  23.  Stress  distribution  across  a  12  cm  thick  steel  plate. 
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Figure  24.  Displacements  across  a  4  cm  thick  steel  plate. 
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Figure  25.  Displacements  across  a  6  cm  thick  steel  plate. 
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Figure  26.  Displacements  across  an  8  cm  thick  steel  plate. 
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Figure  27.  Displacements  across  a  10  cm  thick  steel  plate. 
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Figure  28.  Displacements  across  a  12  cm  thick  steel  plate. 
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Max  Von  Mises  Stress  &  Deflection  vs.  Thickness 


50 

45 

40 

35 

30 

25 

20 

15 

10 

5 

0 


■  Max  Von  Mesis  Stress 


-  Max  Deflection 


Figure  29.  Chart  of  steel’s  performance  at  different  thicknesses. 

Figure  29  shows  the  effect  of  thickness  on  the  maximum  deflection  &  Von  Mises 
stress  for  the  armor  plate.  It  is  observed  that  as  the  thickness  increases  the  maximum 
stress  moves  from  the  center  to  the  edges.  The  maximum  stress  also  decreases  as 
thickness  increases  from  4cm  to  10  cm  followed  by  an  increase  at  12  cm.  In  all  but  the  4 
and  8  cm  thick  plates  the  maximum  Von  Mises  stress  occurs  at  the  top  and  bottom 
constraint,  while  in  the  8  cm  thick  plate  the  maximum  Von  Mises  stress  occurs  on  the  left 
and  right  side.  In  the  4  cm  thick  plate  the  maximum  Von  Mises  stress  occurs  between  the 
center  and  the  top  and  bottom  constraints.  The  discrepancies  in  the  location  of  the 
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Max  Deflection  [mm] 


maximum  Von  Mises  stress  may  be  attributed  to  the  short  simulation  run  time  and/or  the 
whipping  of  the  armor  plate  as  a  result  of  explosive  shock.  The  increase  in  the  maximum 
stress  at  12  cm  of  thickness  may  indicate  the  limitations  of  2D  plate  analysis,  where  the 
thickness  can  violate  the  assumptions  on  the  lateral  normal  stress  and  transverse  shear 
stresses. 

B.  COMPOSITES  WITH  INTER-LAMINAR  SHEAR 

This  section  presents  analysis  results  for  laminated  composite  design  using  TSAI- 
WU  criteria.  Figure  30  and  32  present  the  stresses  and  strains  for  the  mechanical  loads 
only. 
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Figure  30.  Tsai-Wu  criteria  stresses  in  high-modulus  graphite. 
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Figure  3 1 .  Tsai-Wu  criteria  strains  in  high-modulus  graphite. 


When  Figure  31  is  compared  to  Figures  11  and  12,  high-modulus  and  high- 
strength  graphite  as  an  ortho  tropic  material,  the  resulting  stress  is  reduced  to  about  50% 
of  yield  strength  while  the  graphite  material  had  failed.  Another  significant  change  in  the 
results  is  that  the  stress  and  strain  is  distributed  throughout  the  entire  plate,  taking 
maximum  advantage  of  composites  unique  tailoring  properties  compared  to  metals.  The 
resulting  strains,  seen  in  Figure  31,  are  safely  below  failure  criteria  for  typical 
composites.  This  advantage  may  decrease  with  decreases  in  the  size  of  the  plate. 

C.  THERMAL  LOADING  RESULTS 


The  results  from  the  thermal  analysis  are  presented  in  the  following  section.  The 
thermal  loads  applied  to  the  armor  plate  include  the  ambient  temperature  and  the  initial 
temperature  of  the  EFP. 
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Figure  32.  Temperature  distribution  on  an  aluminum  plate. 
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Figure  33.  Temperature  distribution  on  a  steel  plate. 
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Figure  34.  Von  Mises  stress  distribution  caused  by  initial  temperature  on  aluminum. 
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Figure  35.  Von  Mises  stress  distribution  caused  by  initial  temperature  on  steel. 
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Figure  36.  Displacements  from  initial  temperature  on  aluminum. 
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Figure  37.  Displacements  from  initial  temperature  on  steel. 
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The  results  show  that  the  heat  from  the  EFP  does  not  dissipate  very  much  in  one 
half  second  of  contact.  The  stresses  caused  by  the  thermal  loading  are  of  the  same  order 
of  magnitude  as  the  stresses  caused  by  the  kinetic  loads.  The  resulting  deflections  are 
relatively  small  compared  to  those  caused  by  the  kinetic  loads. 

D.  3-D  THERMAL  LOADING  RESULTS 

The  three  dimensional  plate  was  composed  of  ten  nodded  tetrahedrons  increasing 
the  accuracy  of  the  analysis  by  having  higher  order  approximations  for  the  unknown 
variables  (temperatures  and  displacements). 


Figure  38.  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  0.0005. 
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Figure  39.  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  0.507. 


Figure  40.  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  1 .003. 
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Figure  41 .  Initial  temperature  acting  on  a  3d  aluminum  plate  at  t  =  1 .499. 
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Figure  42.  Heat  flux  acting  on  a  3d  aluminum  plate  at  t  =  0.0005. 
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Figure  43.  Fleat  flux  acting  on  a  3d  aluminum  plate  at  t  =  1 .499. 
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Figure  44.  Temperature  gradient  across  a  3d  aluminum  plate  at  t  =  0.0005. 
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Figure  45 .  Temperature  gradient  across  a  3d  aluminum  plate  at  t  =  1 .499. 


As  shown  in  Figures  38  through  41  the  initial  temperature  drops  from  a  high  of 
1400  K  to  a  low  of  342  K  on  the  surface  in  1.5  seconds  while  raising  the  rest  of  the 
aluminum  armor  plates  temperature  by  11.9  K.  The  EFP’s  point  of  impact  has  to  absorb 
large  amounts  of  energy  as  shown  in  Figures  42  and  43.  Figures  44  and  45  show  how  the 
temperature  gradient  at  the  point  of  impact  increases  rapidly  from  173  k/m  to  1440  k/m. 

E.  3-D  COMBINED  LOADING  OF  ALUMINUM 

A  course  mesh  was  selected  for  the  3-D  analysis  resulting  in  2,000  DOF  (degrees 
of  freedom)  as  compared  to  20,000  DOF.  The  high  fidelity  model  needs  enormous 
computing  resources  for  these  simulations  and  was  not  pursued.  The  following  results 
are  for  a  10  cm  thick  plate  of  aluminum  armor  plate  loaded  simultaneously  with  the 
pressure  wave,  initial  velocity,  and  the  thermal  loading. 
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Figure  46.  Temperature  gradient  across  a  3d  aluminum  plate  at  t  =  0.5  sec. 
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Figure  47.  Maximum  von  Mises  stress  on  a  3d  aluminum  plate  at  t  =  0.5  sec. 
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Figure  48.  Maximum  displacements  on  a  3d  aluminum  plate  at  t  =  0.5  sec. 


As  shown  in  Figures  47  and  48  the  maximum  stresses  and  displacements  are 
similar  to  the  2D  analysis  shown  in  figures  7  and  13.  In  both  cases  the  metal  will  yield 
and  require  an  elastic-plastic  analysis  to  determine  any  further  failure  information.  The 
3D  analysis  does  show  a  more  realistic  model  of  the  distortions  caused  by  the  impact  of 
the  EFP. 
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Figure  49.  Maximum  von  Mises  stress  on  a  3d  steel  plate  at  t  =  0.5  sec. 
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Figure  50.  Maximum  displacements  on  a  3d  steel  plate  at  t  =  0.5  sec. 


The  maximum  Von  Mises  stress  in  the  two  dimensional  model  is  an  order  of 
magnitude  larger  than  the  maximum  stress  in  the  three  dimensional  model  as  shown  in 
Figures  8  and  49.  The  corresponding  maximum  deflections  are  24  mm  and  29  mm 
respectively  as  shown  in  Figures  14  and  50 
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VI.  CONCLUSION 


A  methodology  was  developed  and  demonstrated  to  model  and  simulate 
survivable  armor  plates  subjected  to  specified  requirements.  The  design  constraints 
included  initial  temperature  flux  doe  to  EFP,  initial  velocity  due  to  EFP,  and  a 
mathematical  representation  of  blast  loading. 

The  parametric  design  study  included  variations  of  different  materials,  different 
thicknesses,  and  configurations  that  included  composite  construction. 


The  studies  indicate  that  modeling  the  armor  plate  with  2-D  plate  elements  are 
inadequate  due  to  thickness  effects  and  3-D  elements  capture  the  behavior  of  armor.  This 
can  be  evidenced  in  Table  1  and  Figure  29. 


Max  Stress 

Thickness  Von  Mises  Max  Deflection  Weight 
[cm]  [MPa]  [mm]  [kg] 

4 

3250 

45.5 

1,886 

6 

2090 

25.6 

2,830 

8 

1700 

23.2 

3,773 

10 

1690 

24.1 

4,716 

12 

2250 

11.3 

5,659 

Table  1.  Summary  of  key  data  for  steel  of  various  thickness  x  3  m  x  2  m. 

The  high-strength  graphite  composite  perfonned  the  best  followed  by  the  high- 
modulus  graphite  composite  when  the  weight  is  considered  as  show  in  Table  2.  The 
analyzing  software  however  did  not  consider  the  thermal  stresses  in  composites.  A 
recommendation  would  be  to  use  a  sandwich  construct  of  aluminum/high-strength 


51 


composite  with  aluminum  sheet  exposed  to  the  high  temperature/velocity  impact,  this 
approach  would  take  advantage  of  both  the  material’s  properties  for  the  armor  design. 


Density  Weight 
[kg/m3]  [kg] 

Steel 

7,860 

4,716 

Aluminum 

2,800 

1,680 

High-Mod  Graphite 

1,200 

720 

High-Str  Graphite 

1,400 

840 

Table  2,  Weight  of3mx2mx0.ini  armor  plate  made  of  listed  materials. 

Follow  up  analysis  with  this  EFP  model  could  incorporate  ballistic  ceramic 
materials.  The  methodology  established  could  be  used  to  incorporate  candidate  materials 
and  field  relevant  blast  loads  and  initial  velocity  estimates  to  design  the  armor  to  survive 
and  be  less  vulnerable  to  the  IED  threat. 

Any  armor  successful  at  defending  against  the  insurgents  IED  attacks  will  result 
in  the  insurgents  using  bigger  IEDs.  Even  the  composite  materials  used  here  were  heavy 
and  may  prove  impractical  on  military  vehicles.  In  an  overall  survivability  approach  it 
may  prove  more  cost  effective  to  reduce  the  susceptibility  of  military  vehicles,  with  better 
sensors  and  countermeasures,  vice  building  an  oversized  vehicle  with  decreased 
vulnerability. 
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APPENDIX 


A.  ARMOR  PLATE  CONSTRUCTION 

This  section  will  discuss  how  to  create  the  armor  plate  in  PATRAN  starting  with 
opening  a  new  file,  and  ending  with  applying  the  boundary  conditions. 

1.  Our  First  Step  is  to  Create  a  New  Data  Base 

From  the  top  menu  choose  File,  and  in  the  resulting  pull  down  menu  choose  New. 
A  side  menu  called  New  Database  pops  up: 

•  Turn  on  (check)  Modify  Preferences. 

•  Under  Look  In  select  the  computer  directory  you  want  to  work  in. 

•  Under  File  Name  enter  the  name  you  want  to  give  the  file. 

2.  Next  Set  the  Analysis  Preference 

A  New  Model  Preferences  will  appear  as  a  side  menu  on  the  right  hand  side: 

•  Under  Tolerance  choose  based  on  model. 

•  Set  Dimensions  to  10.0. 

•  Under  Analysis  Code  choose  MSC.Nastran 

•  Under  Analysis  Type  choose  Structural. 

3.  The  Geometry  of  the  Beam  will  Beam  Determined  Next 

From  the  top  menu  choose  the  Geometry  button.  A  side  menu  called  geometry 
will  open  on  the  right  hand  side: 

•  Under  Action  select  create. 

•  Under  Object  select  solid. 
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•  Under  Method  select  XYZ 

•  In  the  Vector  Coordinates  List  box  enter  the  dimension  of  the  annor 
plate.  For  example  <3,  2,  0. 1>. 

•  Uncheck  Auto  Execute. 

•  In  the  Origin  Coordinates  List  enter  the  point  you  want  the  armor  plate 
to  originate  at.  Only  the  z  coordinate  should  change  base  the  how 
layers  the  plate  is  made  of. 

•  Click  Apply 

•  Repeat  the  geometry  steps  as  needed  to  create  multiple  layers. 
Adjusting  the  Origin  Coordinates  List  as  needed. 

4.  Specify  the  Boundary  Conditions  Next 

Lrom  the  top  menu  choose  the  Loads/BCs  button.  A  side  menu  called 
Load/Boundary  Conditions  will  appear  on  the  right  hand  side: 

•  Under  Action  select  create. 

•  Under  Object  select  displacement. 

•  Under  Type  select  nodal. 

•  Set  Current  Load  Case  to  default. 

•  In  the  New  Set  Name  box  enter  the  name  you  want  to  give  the 
boundary  condition. 

•  Click  the  Input  Data  button,  now  enter  <0,0, 0>  in  the  Translations  box 
and  in  the  Rotations  box,  and  click  ok. 

•  Click  the  Select  Application  Region  button,  on  the  mouse  pointer 
select  menu  depress  the  Curve  or  Edge  button.  Now  select  the  edges 
to  be  clamped  so  there  label  appears  in  Select  Geometry  Entities  box 
and  click  add.  When  finished  adding  click  ok. 

•  Click  Apply. 
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5. 


The  Finite  Element  Mesh  is  Specified  Next 


From  the  top  menu  choose  the  Elements  button.  A  side  menu  will  appear  on  the 
right  had  side. 

•  Under  Action  select  create. 

•  Under  Object  select  mesh. 

•  Under  Type  select  solid. 

•  Under  Elem  Shape  select  tet. 

•  Under  Mesher  select  IsoMesh. 

•  Under  T opology  select  Tet  1 0 . 

•  In  the  Solid  List  box  enter  all  solids  by  highlighting  them  with  the 
mouse. 

•  In  the  Global  Edge  Length  Box  make  sure  Automatic  Calculation  is 
checked. 

•  Click  Apply 

If  using  multiple  layers  perform  the  following. 

•  Under  Action  select  equivalence. 

•  Under  Object  select  all. 

•  Under  Method  select  tolerance  cube. 

•  In  the  Equivalencing  Tolerance  box  enter  0.005. 

•  Click  Apply 


55 


B.  ARMOR  PLATE  LOAD  CASE  CONSTRUCTION 

1.  The  Load  Case  is  Created  Next 

From  the  top  menu  choose  Load  Cases.  A  side  menu  will  appear  on  the  right 
hand  side. 

•  Under  Action  select  create. 

•  In  the  Load  Case  Name  box  enter  the  name  of  the  load  case  for 
example  “LC1”. 

•  Click  apply. 

2.  The  Loads  are  Specified  Next 

From  the  top  menu  choose  the  Loads/BCs  button.  A  side  menu  called 
Load/Boundary  Conditions  will  appear  on  the  right  hand  side: 

•  Under  Action  select  create. 

•  Under  Object  select  pressure. 

•  Under  Type  select  element  uniform. 

•  In  the  Current  Load  Case  box  click  the  button  and  select  the  load  case 
you  created. 

•  In  the  New  Set  Name  box  enter  the  name  of  the  pressure  force. 

•  Click  the  Input  Data  button  and  enter, 

i.  Translations  <0,0,45000>. 

ii.  Rotations  <0,0,0>. 

iii.  Click  ok. 

•  Click  Select  Application  Region, 
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i.  In  the  Application  Region  box  enter  the  name  of  the  surface  the 
pressure  is  being  applied  to  by  highlighting  it  in  the  work 
window  with  the  mouse  pointer  set  to  Surface  or  Face. 

ii.  Click  add,  the  name  of  the  surface  will  appear  in  the 
Application  Region  box. 

iii.  Click  ok. 


•  Click  apply. 

The  initial  velocity  is  created  next. 

•  Under  Action  select  create. 

•  Under  Object  select  initial  velocity. 

•  Under  Type  select  nodal. 

•  In  the  Current  Load  Case  box  click  the  button  and  select  the  load  case 
you  created. 

•  In  the  New  Set  Name  box  enter  the  name  of  the  initial  velocity. 

•  Click  the  Input  Data  button  and  enter, 

i.  Translations  Veloc  <0,0,-2000>. 

ii.  Rot  Veloc  <0,0, 0>. 

iii.  Click  ok. 

•  Click  Select  Application  Region, 

i.  In  the  Application  Region  box  enter  the  name  of  the  node  the 
initial  velocity  is  being  applied  to  by  highlighting  it  in  the  work 
window  with  the  mouse  pointer  set  to  Point  or  Vertex. 

ii.  Click  add,  the  name  of  the  node  will  appear  in  the  Application 
Region  box. 

iii.  Click  ok. 

•  Click  apply. 


57 


C.  CREATING  ARMOR  MATERIAL 

1.  The  Materials  are  Specified  Next 

From  the  top  menu  choose  the  Materials  button.  A  side  menu  called  Materials 
will  appear  on  the  right  hand  side: 

•  Under  Action  select  create. 

•  Under  Object  select  isotropic  for  metals  and  2d  ortho  tropic  for 
composites. 

•  Under  Method  select  Manual  input. 

•  In  the  Material  Name  box  enter  the  name  of  the  material  (i.  e.  Steel). 

•  Click  the  Input  Properties  button  for, 

i.  Under  Constitutive  Model  select  Linear  Elastic.  For  metals 
enter  the  values  of  Elastic  Modulus,  Poisson  Ratio,  Shear 
Modulus,  and  Density.  For  composites  enter  the  values  of 
Elastic  Modulus  11,  Elastic  Modulus  22,  Poisson  Ration  12, 
Shear  Modulus  12,  Sheer  Modulus  23,  and  Density. 

ii.  Under  Constitutive  Model  select  Failure.  For  metals  enter  the 
values  of  Tension  Stress  Limit,  Compression  Stress  Limit,  and 
Shear  Stress  Limit.  For  composites  select  Failure  1  and  enter 
the  values  of  Tension  Stress  Limit  1 1,  Tension  Stress  Limit  22, 
Compress  Stress  Limit  11,  Compress  Stress  Limit  22,  Shear 
Stress  Limit. 

iii.  Click  ok. 

•  Click  apply. 
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Tension 

480 

480 

impression 

480 

480 

Shear 

480 

290 

Modulus  of  Elasticty  [Gpa] 

El 

200 

72 

290 

128 

E2 

6.9 

11 

Modulus  of  Rigidity  [Gpa] 

G 

79 

28 

4.8 

4.5 

Piosions's  Ratio 

Nu 

0.30 

0.30 

0.25 

0.25 

Density  [kg/m3] 

Ro 

7860 

2800 

1200 

1400 

Yield  Strength  [Mpa] 

Tension 

345 

410 

Xt 

621 

1170 

Xc 

621 

1120 

Yt 

13.8 

41 

Yc 

193 

172 

Shear 

210 

220 

S 

27.6 

48 

Coefficient  of  Thermal  Expansion  [10A- 

-6  C] 

alphal 

11.7 

23 

-0.32 

0.14 

alpha2 

9.17 

8.44 

Specific  Heat  [J/kg  *  K] 

Cp 

446 

875 

935 

935 

Thermal  Conductivity  [W/m 

*  K] 

k  parallel 

51.9 

177 

11.1 

11.1 

k  perpendicular 

0.87 

0.87 

Table  3.  Properties  of  various  tested  materials. 
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Creating  a  Composite 

While,  still  working  under  the  materials  menu: 

•  Under  action  select  create. 

•  Under  object  select  composite,  a  Laminated  Composites  window  will 
appear. 

•  In  the  insert  box  enter  the  number  of  layers  to  be  created. 

•  For  each  row  click  the  box  under  the  material  heading  and  select  a 
material  from  the  materials  list  showing  on  the  materials  right  hand  menu. 

•  Next  click  in  the  box  under  the  thickness  heading  enter  the  thickness  in  the 
in  the  input  data  box. 

•  Do  the  same  for  the  box  under  the  orientation  heading. 

•  When  done  entering  laminate  layers  click  the  click  the  show  laminate 
properties  button,  and  the  properties  of  the  laminate  will  be  displayed  in 
matrix  form. 

•  To  close  the  laminate  select  another  option  under  object  or  choose  another 
top  menu  option. 

2.  The  Properties  are  Assigned  Next 

From  the  top  menu  choose  the  Properties  button.  A  side  menu  called  Properties 
will  appear  on  the  right  hand  side: 

•  Under  action  select  create. 

•  Under  object  select  2D. 

•  Under  type  select  shell. 

•  In  the  Property  Set  Name  box  enter  a  name. 

•  Click  the  Input  Properties  button,  an  input  properties  window  will  apper. 
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i.  Click  the  cross  hatch  icon  next  to  Mat  Prop  Name  and  select  a 
material  from  the  list  that  appears. 

ii.  For  2D  objects  enter  the  thickness  in  the  thickness  box 

iii.  Click  ok. 

•  Click  the  Select  Application  Region  button,  a  new  menu  will  open  on  the 
right. 

i.  Click  in  the  Select  Members  box,  then  go  to  the  work  space  and 
select  the  whole  annor  plate. 

ii.  Click  the  add  button,  the  name  of  the  whole  plate  will  appear  in  the 
Application  Region  box. 

iii.  Click  ok. 

•  Click  Apply 

D.  KINEMATIC  LOAD  ANALYSIS 
1.  Analysis  is  Performed  Next 

From  the  top  menu  choose  the  Analysis  button.  A  side  menu  called  Analysis  will 
appear  on  the  right  hand  side: 

•  Under  action  select  analyze. 

•  Under  object  select  entire  model. 

•  Under  method  select  analysis  deck. 

•  In  the  job  name  box  enter  a  job  name. 

•  Click  the  Translation  Parameters  button. 

i.  Select  OP2,  and  click  ok. 

•  Click  the  Solution  Type  button. 


61 


i.  Select  Transient  Response  for  the  solution  type. 

ii.  Under  Fonnulation  select  direct. 

iii.  Click  ok. 

•  Click  the  subcases  button. 

i.  In  the  Available  Subcases  box  highlight  the  subcase  to  be 
analyzed,  the  name  will  appear  in  the  Subcase  Name  box. 

ii.  Highlight  the  load  case  to  be  analyzed  in  the  Available  Load  Cases 
box. 

iii.  Click  the  Subcase  Parameters  button,  then  click  the  Define  Time 
Steps  button. 

iv.  Enter  the  number  of  time  steps  and  the  value  of  delta  t,  then  click 
ok. 

v.  Click  ok  on  the  Subcase  Parameters  box. 

vi.  Click  the  Out  Put  Requests  box,  and  select  the  desired  outputs 
form  the  menu,  and  click  ok. 

vii.  Click  Apply  in  the  Subcases  window,  then  cancel  to  close  the 
window. 

•  Click  the  subcase  select  button. 

i.  Highlight  the  desired  load  case  in  Subcases  For  Solution  Sequence 
box. 

ii.  Remove  any  extra  load  cases  from  the  Subcases  Selected  box. 

iii.  Clock  ok. 

•  Click  Apply,  PATRAN  will  create  the  .bdf  file  in  temp  directory. 

•  Next  open  Nastran,  select  the  .bdf  file  and  click  open. 

•  In  the  option  keywords  box  enter,  “news=no,”  and  ,  “scr=yes”. 

•  Click  Run. 
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2. 


Accesses  the  Results 


Open  a  database  file.  Then  click  the  Analysis  button. 

•  Under  action  select  access  results. 

•  Under  object  select  read  output2. 

•  Under  method  select  both. 

•  Click  apply. 

•  Then  from  the  top  menu  click  Results. 

•  View  the  results  creating  images  or  videos  of  desired  data. 

E.  THERMAL  LOAD  ANALYSIS 

From  the  top  menu  choose  the  Loads/BCs  button.  A  side  menu  called 
Load/Boundary  Conditions  will  appear  on  the  right  hand  side: 

•  Under  Action  select  create. 

•  Under  Object  select  initial  temperature. 

•  Under  Type  select  element  uniform. 

F.  THERMAL  STRESS  ANALYSIS 

To  perform  a  stress  analysis  you  will  need  open  a  new  data  base,  and  click  the 
Analysis  button  on  the  top  menu. 

1.  Accesses  Thermal  the  Results 

•  Under  action  select  access  results. 

•  Under  object  select  read  output2. 

•  Under  method  select  both. 


63 


•  Click  Results  on  the  top  menu,  and  select  any  result  to  view,  you  must  do 
this  for  a  later  state. 

2.  Create  a  Spatial  FEM  Based  on  the  Temperature  Profile 

•  Under  action  select  create. 

•  Under  object  select  spatial. 

•  Under  method  select  FEM. 

•  In  the  Field  Name  box  enter  a  name. 

•  In  FEM  file  definition  select  continuous. 

•  In  Mesh/Results  group  filter  select  Scalar. 

•  In  select  group  select  default  group,  if  you  did  choose  a  set  of  result  this 
will  be  blank. 

•  Click  apply 

3.  Change  the  Analysis  Type 

•  Select  Preferences/ Analysis  for  the  menu  above  the  top  menu. 

•  Select  Analysis  type,  choose  structural,  and  click  ok. 

•  Add  or  change  the  materials  as  before  if  needed. 

•  Add  or  change  the  properties  as  before  if  needed. 

•  Create  a  new  load  case  as  before,  the  constraint  forces  will  need  to  be 
added  or  the  temperature  distribution  can  be  added  to  an  existing  load 
case. 

4.  Define  a  Temperature  Load 

•  Under  action  select  create. 
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•  Under  object  select  temperature. 

•  Under  type  select  Nodal. 

•  In  the  New  Set  Name  box  enter  the  name  of  the  temperature  load. 

•  Click  the  input  data  button,  and  select  the  name  to  the  temperature  from 
the  thermal  analysis,  and  click  ok. 

•  Click  the  select  application  region  button,  choose  the  whole  plate  of 
annor,  and  click  ok. 

•  Click  Apply 

•  Perform  the  analysis  as  before. 
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